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Abstract 

The neutron flux induced by radioactivity at the Gran Sasso underground labora- 
tory is revisited. We have performed calculations and Monte Carlo simulations; the 
CN \ results offer an independent check to the available experimental data reported by 

different authors, which vary rather widely. This study gives detailed information 
on the expected spectrum and on the variability of the neutron flux due to possible 
variations of the water content of the environment. 
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1 Introduction 



One of the important parameters characterizing an underground environment 
is the neutron flux. Recently, the role of neutrons as a source of background 
has received more attention. Estimates of the neutron flux from surrounding 
materials are needed for rare event searches in underground sites. Highly sensi- 
tive underground neutrino experiments require precise knowledge of neutron- 
induced background in the detector. In direct searches for WIMP dark matter 
neutrons are a particularly important background, because they interact with 
the detector in the same way WIMPs do. 

The dominant sources of neutrons at large depth underground, where the 
cosmic rays are reduced significantly, are the (a,n) reactions on light elements 
(e.g. Li, F, Na, etc.) and the spontaneous fission, mainly of 238 U. In general, 
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at a depth of 3000 — 4000 m w.e., where many rare event experiments are 
located, the flux of neutrons from activity of the experimental environment 
is two to three orders of magnitudes higher than the flux of neutrons from 
cosmic ray muons. 

Due to its very low intensity, the neutron flux underground is not easy to 
measure. At the Gran Sasso underground laboratory the neutron flux has 
been measured by several groups, employing different methods of detection 
[1,2,3,4,5,6]. To get the neutron flux from the measurements, Monte Carlo 
simulations are usually performed, because the parameters required to handle 
the experimental data (detector efficiency, energy spectrum of the incoming 
neutrons, etc.) are usually unknown or difficult to determine experimentally. 
Table 1 summarizes the results of the neutron flux measurements at the Gran 
Sasso laboratory. One measurement was performed in hall C by the ICARUS 
collaboration [2], while all others were performed in hall A of the laboratory. 
Most of the measurement results are given as the integral flux in a large energy 
bin. Although different energy binning was used in these measurements, one 
can see that the results vary rather widely. Only two measurements [2,3] give 
some information on the spectral shape. However, for some purposes a more 
detailed spectrum is needed. 

The total flux of [2] is basically consistent with the result of [3] for E > 1 MeV, 
while the energy spectrum as a whole is significantly harder. The results from 
[2] show important distortions with respect to the spectrum that was pre- 
viously generally assumed in literature, i.e. a spectrum produced mainly by 
spontaneous fission. The difference is seen at high energies. To search for a 
possible explanation the authors of [2] considered the production of neutrons 
via (a,n) reactions inside the rock and computed the 'thin target' spectrum 
of neutrons produced by the reaction 17 O(a:,n) 20 Ne only. The (a,n) and spon- 
taneous fission spectra were then used as input for the MCNP (Monte Carlo 
N-Particles) code to obtain the effective neutron flux in hall C. They found 
that their measurements and simulations are in very good agreement and that 
(a,n) reactions are the main source of the high energy neutron flux in the lab- 
oratory. However, according to [7] the uncertainties quoted in [2] have been 
underestimated. They should be multiplied by \/3 to get the proper uncer- 
tainties, due to the fact that the background run in this experiment is about 
three times shorter than the data run. Then all data points shown in [2] would 
be insignificant and the agreement between the thin target calculation and the 
measurement found here would be merely accidental. 

Based on Monte Carlo simulations, this work is devoted to reinvestigate the 
neutron flux from (a,n) and fission reactions in the rock and concrete of the 
Gran Sasso laboratory. The results offer an independent check to the present 
experimental data and give detailed information on the shape of the spectrum. 
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2 Calculation of Neutron Production by («,n) and Fission Reac- 
tions in the Rock and Concrete 

2. 1 Composition and Activity of the Rock and Concrete 

Gran Sasso rock consists mainly of CaC0 3 and MgC0 3 , with a density of 
2.71±0.05 g/cm 3 [8]. The weight percentage of the elements is given in Table 2. 
Due to the presence of a certain type of rock, called "roccia marnosa nera" , 
the contaminations of 238 U and 232 Th in hall A rock are about ten and thirty 
times higher respectively than those of hall C [4] as shown in Table 3. 

Since there are no data on the chemical composition of Gran Sasso concrete 
available in literature, several samples were taken from different places in the 
laboratory. The concrete samples were then analyzed at the laboratory of the 
Lehrstuhl fur Baustoffkunde und Werkstoffpriifung at the department of civil 
engineering of the Technische Universitat Munchen. No significant variations 
were found in the chemical composition of the samples, which leads to the 
conclusion that all halls in Gran Sasso are layered with the same type of 
concrete. The typical water content in the concrete is 12%, with a possible 
variation of 4% at most (in most cases the variation is smaller). The weight 
percentage of elements in concrete with 8% water content (hereafter "dry 
concrete") is shown in Table 4. The 238 U and 232 Th contaminations are 1.05 ± 
0.12 ppm and 0.656 ± 0.028 ppm [9] respectively. The density is between 2.3 
and 2.5 g/cm 3 , depending on the assumed water content. 

2.2 Neutron Production by Spontaneous Fission 

There are mainly three nuclides in nature that undergo spontaneous fission: 
238 U, 235 U and 232 Th. Only neutrons produced by spontaneous fission of 238 U 
are considered here, because of the long fission half life of the other two nuclides 
compared to that of 238 U. The spectrum of the emitted neutrons follows the 
Watt spectrum: 

N(E) = Cexp (-£/a)sinh {bE) 1/2 (1) 

In this work the Watt spectrum parameters of the Los Alamos model results 
[10] have been used, where a = 0.7124 MeV and b = 5.6405 MeV" 1 . The rate 
of spontaneous fission of 238 U is 0.218/year/g of rock (concrete) for lppm of 
238 U and the average number of neutrons emitted per fission event is 2.4±0.2 
[11]. This gives 0.52 neutrons/year /g of rock (concrete) /ppm 238 U. This num- 
ber multiplied by the 238 U activities in the rock/concrete would give 3.54, 
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0.22 and 0.34 neutrons/year /g in the rock of hall A, B and C respectively 
and 0.55 neutrons/year /g in the concrete. Figure 1 shows the neutron energy 
distribution due to spontaneous fission of 238 U. 



2.3 Neutron Production by (a,n) Reactions 



Uranium, thorium, and their daughter products decay by emitting a and (5 
particles. In the rock (concrete) a-particles can interact especially with light 
elements and produce neutrons through (a,n) reactions. 

The yield of neutrons per a-particle for an individual element depends on the 
(cc,n) interaction cross section (which is energy dependent), and on the energy 
loss of a-particles in a medium made of that element. In this work the thick 
target yield of (a,n) reactions was used instead of the thin target yield used 
in [2]. The thick target yield of the (a,n) reaction for an individual element j 
in which the a-particle has a range R can be written as [12]: 

R 

Yj = Jnja^dx (2) 
o 



where rtj is the number of atoms per unit volume of element j, and o~j is the mi- 
croscopic (a,n) reaction cross section for an a-particle energy E. Transforming 
the right side of Eq. 2 into an integral over energy gives: 



Y >-l -{dE/dx) dE 



o 

E ■ 

E ' n i a j{E)NA 
njAjSfiE) 



dE 



N *f"'^.iE (3) 



Aj J SJ-(E) 

where E^ is the initial a-energy, Na is Avogadro's number, Aj is the atomic 
mass and <7j and S™ are the (a,n) cross section and the mass stopping power 
respectively, which are energy dependent. 

Neutron yields from individual elements can be used to calculate the total yield 
in a chemical compound or mixture [12,13,14]. The following assumptions are 
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usually made in such a calculation: 

(i) the compound is a homogeneous mixture of its constituent elements 

(ii) Bragg's law of additivity for stopping power holds for the compound 

(iii) the ratio of an element's stopping power to the total stopping power of 
the compound is independent of the a-particle energy. 

The validity of (iii) has been discussed by several authors. Although for each 
element the mass stopping power S™ decreases with energy, Feige [12] found 
that the mass stopping power ratio for any pair of elements does not change 
by more than ±4% between 5.3 MeV and 8.8 MeV. Heaton et al. [15] showed 
that above 3 MeV this approximation introduces an uncertainty of less than 
5% in the neutron yield. 

Under those assumptions the neutron yield of element j in the compound or 
mixture with initial a-particle energy Ei can be written as: 



Ei 

- J2jPjS'j 



_ r njaj 



dE 



P 3 Sf{E ) (n 3 a,{E) 



E jPj S™(e )J Pj S) 



ZjMjS^Eo) 



Y,(E t ) (4) 



where Mj is the mass fraction of element j in the mixture, E is a chosen 
reference energy (8 MeV in this work), S™ is the mass stoping power and 
Yj(Ei) is the neutron yield of element j in isolation (see Eq. (3)). Thus the (a,n) 
yield of a compound (mixture) is the sum of the yields of its elements weighted 
by the relative contributions of the elements to the total stopping power of 
the compound. The use of 8 MeV mass stopping power was selected because 
the overwhelming contribution to the neutron yield comes from high energy 
a-particles and the relative stopping power of elements are nearly independent 
of the a-energy in this energy region [16] . The mass stopping powers at 8 MeV 
of elements used in this work are as those used in [15]. In this work the data 
on neutron yield of elements in isolation, in units of neutron/a, is taken from 
a compilation by Heaton et al. [16,17]. 

Each a emitter in the 238 U and 232 Th decay chains emits cc's at a certain energy, 
which was used as the initial energy in the neutron yield calculation. The 
neutron yield of each element with certain initial energy was then multiplied 
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by the branching ratio, the number of ct's emitted by each emitter per unit 
time, and the concentration of 238 U and 232 Th in the rock (concrete). It is 
assumed in this work, that 238 U and 232 Th are in secular equilibrium with 
their daughter products. The total neutron production rate for each element 
was calculated by summing up the contribution of all alpha emitters. 

In Table 5 and Table 6 neutrons produced per unit mass per year by (a,n) 
reactions with the elements of the rock in hall A, hall C and also in dry concrete 
are presented. It is seen that fission (see discussion in the previous section) and 
(a,n) reactions contribute more or less equally to the total production rates 
both in the rock and in the concrete. Half of the total (a,n) neutron production 
in the rock comes from interactions of a particles with magnesium, which 
comprises only less than 6% of the weight percentage of the rock, whereas 
oxygen with almost 50% weight percentage contributes to only about 20% of 
the production rate. Due to the higher activity of the hall A rock the (a,n) 
neutron production in the rock of this hall is more than ten times higher than 
in the hall C rock. In the concrete Na, Al and Mg contribute significantly in 
spite of their minor weight percentages. The production rate per unit mass 
in the wet concrete (16% water content) is slightly smaller than in the dry 
concrete, but this is merely due to the difference in the densities. The volume 
of concrete remains, while the mass changes with the water content. Given 
per unit volume the production rates in dry and wet concrete are the same. 

The energy of the emitted neutron is dependent on the a energy, the reac- 
tion energy Q, and the neutron emission angle. It was calculated under the 
following assumptions: 

(i) the interaction take place at the initial a energy 

(ii) the neutron is emitted at 90° 

(iii) the residual nucleus is produced in its ground state. 

Under these assumptions, the neutron energy can be determined by using the 
simple Eq. (5): 



Here M is the mass of the final nucleus, M n and M a are the masses of neutron 
and the a particle respectively, and E a is the initial a energy. 

The threshold energy E t h for an (a,n) reaction is the minimum kinetic energy 
the impinging a particle must have (in the laboratory system) in order to make 
the reaction energetically possible. For endothermic reactions, the threshold 
energy is: 



MQ + E a (M - Mg) 
(M n + M) 



(5) 



E th = -[(M n + M a )/M 1 ]Q 



(6) 
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where Mi and M a are the masses of the target nucleus and of a particle re- 
spectively. The threshold energy is zero for exothermic reactions, i.e. reactions 
with positive Q. 

The highest energy among naturally emitted relevant a particles is 8.79 MeV, 
which comes from the decay of 212 Po. Hence, for some elements in the rock/ con- 
crete there are isotopes that can not participate in the (a,n) reactions because 
of their high E t h (e.g. 16 0, 28 Si, 24 Mg, and 40 Ca). For each element calculations 
of neutron energies were therefore done for all isotopes that are not closed for 
(a,n) interactions. Then, neutron mean energies of elements were calculated 
according to the relative abundances of the "open" isotopes for different a en- 
ergies. Finally, the yields of all elements were summed up in 0.5 MeV energy 
bins to get the energy spectra of neutrons from (ct,n) reactions as shown in 
Figure 2 for hall A rock, hall C rock, and dry concrete. While neutrons below 
4 MeV are mainly produced by spontaneous fission, (ct,n) reaction is the main 
contributor in the production of neutrons with higher energy. Neutrons with 
energy above 6 MeV are contributions of magnesium and carbon. 



3 Flux of (a,n) and Fission Neutrons at LNGS 

To get the flux of low energy neutrons inside the halls, the Monte Carlo code 
MCNP4B (Monte Carlo N-Particles version 4B) from Los Alamos [18] was used 
to transport neutrons produced by fission and (a,n) reactions as calculated 
above through the rock and concrete and scatter them inside the halls. 

Table 7 summarizes the typical depths of neutrons (the depths where ~ 63% 
of neutrons come from, the 1/e length) produced in hall A rock and in dry 
concrete, and entering hall A (before scattering inside the hall) with any energy 
and in addition with energy E > 1 MeV. Neutrons with E > 1 MeV come 
mainly from the first 7 cm and 13 cm of concrete and rock respectively. As 
the thickness of the concrete layer at the Gran Sasso laboratory is not less 
than 30 cm (in some places even around 1 m), the bulk of the total flux at the 
laboratory is given by neutrons produced in the concrete. 

In the further simulations performed to get the neutron flux inside the halls, 
the thickness of the concrete layer in the laboratory has been set to 45 cm 
below the floor and 35 cm elsewhere, and 2 m of rock have been taken into 
account. The calculated fluxes of neutrons (after scattering inside the halls) 
are shown in Table 8 for hall A with 8% (dry) and 16% (wet) water content in 
the concrete and for hall C with dry concrete, together with the measurement 
in Hall A by Belli et al. [3] as a comparison. The measured flux between 
1 — 2.5 MeV is taken from [2]. We calculated the errors in this energy bin by 
assuming that the relative errors in 1 keV — 1 MeV and 1 — 2.5 MeV bins are 
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equal. 



Uncertainties in the neutron flux in our calculations can be attributed to the 
uncertainties in the neutron production rate from (a, n) /fission reactions and 
the energies of the emitted neutrons. 

The major sources of uncertainty in the total (a,n) neutron yield of the rock 
and concrete are the uncertainties in the stopping powers, elemental neutron 
yield per a particle and material compositions. We adopt the 5% uncertainty 
in the stopping power data as stated in [16]. The assumption that the ratio of 
an element's stopping power to the total stopping power of the compound is 
independent of the a-particle energy introduces another uncertainty of about 
5% in the neutron yield (see Subsection 2.3). The deviation from Bragg's 
Rule, which was observed in compounds in gaseous phase and hydrocarbons, 
does not produce significant contribution to the uncertainty of the neutron 
yield calculation (see [16] and references therein). An additional 5% uncer- 
tainty coming from the Heaton's compilation of elemental neutron yield per 
a above 3.7 MeV [16] is taken into account. For the determination of the 
material composition we assume an uncertainty of 1-2%, as is common using 
modern analytical instruments. Combining these uncertainties in quadrature, 
the uncertainty of the (a,n) yield from a compound is roughly 9%. The major 
uncertainty in the fission neutron yield of about 8% is due to uncertainty in 
the average number of neutrons emitted per fission [11]. The total number of 
neutrons produced by fission and (a,n) in the rock/concrete at the Gran Sasso 
laboratory depends eventually on the 238 U and 232 Th contamination. Taking 
into account the uncertainties in the activities and the aforementioned uncer- 
tainties in the neutron yield would give an over all uncertainties of about 15% 
and 20% in the total neutron fluxes in hall A and hall C respectively, due to 
uncertainties in the neutron production from (a, n) /fission reactions. 

Jacobs and Liskien [19] reported energy spectra of neutrons produced by a- 
particles with energies of 4, 4.5, 5 and 5.5 MeV in thick targets of light elements 
(all elements important for us are available here, except Na). The difference 
between our calculation of the energies of emitted neutrons and the average 
energies in [19] are generally less than 0.5 MeV. Only for Mg a difference of 1.5 
MeV is found. To estimate the uncertainty in the neutron flux due to our as- 
sumptions in determining the energies of the emitted neutrons (see Subsection 
2.3 and Figure 2) we have repeated the simulations of neutron propagation, 
using the average neutron energies reported in [19] where available. For a- 
energies above 5.5 MeV we reduced our calculated emitted neutron energies 
by 0.5 MeV (1.5 MeV for Mg). For elements which are not available from [19] 
we reduced the energies of emitted neutrons by 0.5 MeV for all a-energies. For 
the spontaneous fission of 238 U no modification is needed since the spectrum 
is well understood. We find that the total neutron fluxes at the hall A and 
C change by about 10% due to these alterations, and also the spectral shape 
changes. Combining this uncertainty with the aforementioned uncertainties in 
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the neutron production, we estimate the over all systematic uncertainty of the 
total flux to be around 20% (slightly higher for hall C). Energy dependent 
uncertainties are reported in Table 8. 



The total flux and the integral flux above 1 MeV in hall A are consistent with 
those of the measurement [3] if the concrete in hall A is dry. The simulated 
flux differs from the measured one in the chosen bins below 1 MeV. In our 
simulations we have assumed that there is nothing inside the hall and neutrons 
can only be scattered by the walls. In reality as in the measurement, neutrons 
coming from the rock/concrete can be scattered by anything inside the hall 
before they eventually come into the experimental setup. Those neutrons are 
moderated, raising the flux in the lowest energy bin. 



Within the estimated uncertainties the total flux in hall C is only slightly less 
than in hall A for the case of dry concrete, although the neutron production 
rate in hall C rock is more than ten times lower than that of hall A; above 
1 MeV the fluxes in the two halls are in agreement. This is due to the concrete, 
which indeed reduces the neutron flux from the rock significantly so that 
neutrons coming into the halls are mainly those produced in the concrete 
layer. 



Table 8 shows that the neutron flux depends on the humidity of the environ- 
ment. The flux in Hall A is lower if the concrete is wet than if it is dry (8% and 
16% water content respectively). As mentioned in the previous section, the 8% 
difference in the water content of concrete does not lead to different neutron 
production rates. The effect seen in the flux here is caused only by moderation. 
Wet concrete moderates neutrons more effectively than dry concrete due its 
higher hydrogen content. The fluxes obtained for dry and wet concrete here 
show the maximum possible variation for the water content of concrete. A 
more realistic variation of the water content of (12 ± 1)% results consequently 
in smaller flux variation. To quantify this effect and to see whether it is a sea- 
sonal phenomenon, it is necessary to monitor the water content of the concrete 
for at least one year. 



Detailed spectra of neutrons in hall A and hall C are shown in Figure 3 for 
neutron energy above 0.5 MeV. Each point shows the integral flux in a 0.5 
MeV energy bin. The contribution of (a,n) makes the spectra in both halls 
differ from the spectrum expected for neutrons produced by fission reactions 
only, as was previously generally assumed, especially at high energies. 
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4 Conclusion 



We have discussed the flux of neutrons induced by radioactivity in the rock 
and concrete surrounding the Gran Sasso laboratory. The flux is dominated by 
neutrons produced in the concrete layer and therefore does not vary much from 
hall to hall. It can be expected that as well for other underground laboratories 
the neutron flux originates mainly from the concrete and not from the rock 
material. A more detailed spectrum compared to that from measurements has 
been obtained. The spectrum differs from the spectrum expected for neutrons 
produced by fission reactions only, especially at high energies, due to the 
contribution of (a,n) neutrons. We also have shown the dependence of the 
neutron flux on the humidity of the concrete. Our results for the case of hall 
A with dry concrete are in good agreement with the experimental data from 
[3]. 
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Table 1 

Neutron flux measurements at the Gran Sasso laboratory reported by different 
authors. In analyzing their experimental data with Monte Carlo simulations, Belli 
et al. [3] have used two different hypothetical spectra: flat, and flat plus a Watt 
fission spectrum. This leads to the upper and lower data sets shown for ref. [3] 
respectively. 



E interval 


Neutron Flux (10 6 cm 2 s 1 ) 


(MeV) 


Ref. [1] 


Ref. [2] 


Ref. [3] 


Ref. [4] 


Ref. [5] 


Ref. [6] 


10- 3 -0.5 














u.o — 1 






0.54±0.01 








1 - 2.5 




n i a 4~n 1 o 
U.l4±U.lz 


^U.oo±U.Uo J 








z.o — O 




0.13±0.04 


0.27±0.14 








3-5 






(0.18±0.04) 






2.56±0.27 


5-10 




0.15±0.04 


0.05±0.01 
(0.04±0.01) 


3.0±0.8 


0.09±0.06 




10-15 


0.78±0.3 


(0.4 ±0.4)- lO" 3 


(0.6 ±0.2) TO- 3 
((0.7±0.2)-10~ 3 ) 








15-25 






(0.5±0.3)-10- 6 
((0.1 ±0.3)-10- 6 ) 









Table 2 

Chemical composition of LNGS rock. 



Element 


C 


O 


Mg 


Al 


Si 


K 


Ca 


% Weight 


11.88 


47.91 


5.58 


1.03 


1.27 


1.03 


30.29 



Table 3 

238 U and 232 Th activities in LNGS rock. 



Hall 


Activit 


tes (ppm) 




238tj 


232^ 


A 


6.80 ±0.67 


2.167 ±0.074 


B 


0.42 ±0.10 


0.062 ±0.020 


C 


0.66 ±0.14 


0.066 ±0.025 



Table 4 

Chemical composition of LNGS dry concrete. 



Element 


H 


C 


O 


Na 


Mg 


Al 


Si 


P 


S 


K 


Ca 


Ti 


Fe 


% Weight 


0.89 


7.99 


48.43 


0.6 


0.85 


0.9 


3.86 


0.04 


0.16 


0.54 


34.06 


0.04 


0.43 
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Table 6 



Neutron yields from (a,n) interactions in 
dry concrete (8% water content). 



Table 5 

Neutron yields from ( 
tions in the rock. 


a, n) interac- 


Element 


Total elemental yield 
(n/y/g concrete) 


Element 


Total elemental yield 




C 


5.24E-2 




(n/y/g rock) 







1.50E-1 




Hall A 


Hall C 




Na 


9.65E-2 


C 


4.60E-1 


4.07E-2 




Me 


6.07E-2 





8.8E-1 


7.90E-2 




Al 


5.35E-2 


Mg 


2.31E+0 


2.04E-1 




Si 


3.16E-2 


Al 


3.50E-1 


3.05E-2 




P 


4.35E-4 


Si 


6.00E-2 


5.21E-3 




S 


3.36E-4 


K 


9.00E-2 


7.60E-3 




K 


8.31E-3 


Ca 


2.40E-1 


2.05E-2 




Ca 


4.95E-2 


Total yield 


4.38E+0 


3.88E-1 




Ti 


3.35E-4 








Fe 


9.53E-4 








Total yield 


5.05E-1 



Table 7 

Typical depths (in cm) of (a,n) and fission neutrons produced in hall A rock and 
dry concrete. "Energy" represents the energy of neutrons entering hall A before 
scattering in the hall. 



Energy 


Hall A rock 


Dry concrete 


Fission 


(a,n) 


Fission 


(a,n) 


all energies 


23 


24 


10.5 


12 


> 1 MeV 


10 


13 


6.5 


6.5 
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Table 8 

Simulated neutron flux in hall A with dry (8% water content) and wet (16% water 
content) concrete, in hall C with dry concrete and the measurements by Belli et 
al. [3] for the two hypothetical spectra mentioned in the caption of Table 1. 



Energy 
(MeV) 


Neutron Flux (10 6 n/cm 2 s) 


Measurement in hall A [3] 


MC Simulations, this work 


Flat 


Flat+Watt Spect. 


Hall A, dry 


Hall A,wet 


Hall C,dry 


0- 5.1(T 8 


1.08±0.02 


1.07±0.05 


0.53±0.36 


0.24±0.15 


0.24±0.17 


5.1CT 8 - 1(T 3 


1.84±0.20 


1.99±0.05 


1.77±0.45 


0.71±0.19 


0.93±0.33 


1(T 3 - 2.5 
(1 - 2.5) 


0.54±0.01 
f0.38i0.01) 1 


0.53±0.08 
f0.38i0.06) 1 ' 2 


1.22±0.32 
f0.35±0.12) 


0.57±0.16 
f0.18±0.06) 


0.91±0.32 
f0.27±0.12) 


2.5 - 5 


0.27±0.14 


0.18±0.04 


0.18±0.05 


0.12±0.04 


0.15±0.05 


5-10 


0.05±0.01 


0.04±0.01 


0.05±0.02 


0.03±0.02 


0.03±0.01 


Total Flux 
Flux(E> 1 MeV) 


3.78±0.25 
0.70±0.14 


3.81±0.11 
0.60±0.07 


3.75±0.67 
0.58±0.13 


1.67±0.29 
0.33±0.07 


2.26±0.49 
0.45±0.13 



1 Taken from [2]. Uncertainties are assumed to be equal below and above 1 MeV. 
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Fig. 1. Energy spectrum of neutrons from spontaneous fission of 38 U. 
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Fig. 2. Energy spectra of emitted neutrons from (a,n) reactions as used in the 
simulations. 
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Fig. 3. Neutron flux at the Gran Sasso laboratory, •: hall A, dry concrete, x: hall 
A, wet concrete, 0: hall A, dry concrete, fission reactions only and O: hall C, dry 
concrete. Each point shows the integral flux in a 0.5 MeV energy bin. 
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